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We study the effect of random magnetic fields in the spin-flavor precession of solar neutrinos in 
a three generation context, when a non-vanishing transition magnetic moment is assumed. While 
this kind of precession is strongly constrained when the magnetic moment involves the first family, 
such constraints do not apply if we suppose a transition magnetic moment between the second and 
third families. In this scenario we can have a large non-electron anti-neutrino flux arriving on Earth, 
which can lead to some interesting phenomenological consequences, as, for instance, the suppression 
of day-night asymmetry. We have analyzed the high energy solar neutrino data and the KamLAND 
experiment to constrain the solar mixing angle, tan#o, and solar mass difference, Am%, and we 
■ have found a larger shift of allowed values. 
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O ! I. INTRODUCTION 

Recent results of KamLAND experiment Q confirmed the Large Mixing Angle (LMA) realization of the MSW 
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phenomenon as the explanation to the solar neutrino anomaly l^liUjISllflljjjIalSl- Furthermore, KamLAND results 



rule out several other possible solutions to the solar neutrino problem based on exotic phenomena [lOj], like as the 
^ I resonant spin-flip conversion QIlllQQQQHQ 
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induced by a non-vanishing neutrino magnetic moment 

interacting with solar magnetic fields. 



CIh! Nevertheless exotic phenomena can generate sub-leading effects which are still allowed by present solar neutrino 



data. Such effects can add new features to this picture, in particular, changing the determination of the neutrino 
oscillation parameters. Examples of these sub- leading effects were analyzed in |l8j, where random fluctuations of solar 
matter were considered, or in [ltl Lfj. where non-standard neutrino interactions induced a different determination 
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of the oscillation parameters necessary for a solution to the solar neutrino anomaly. Here we study another possible 
sub- leading effect: the consequences of neutrino interaction with random solar magnetic fields through a non- vanishing 
magnetic moment. 

The random magnetic scenario was studied in a different context l^. I23I I24I I25I l^rl but always assuming 
a magnetic moment linking the electron neutrino with the muon and tau anti-neutrino families. In this framework 
electron anti-neutrinos are produced as a consequence of the spin-flavor precession, due the large mixing angles Oq 
and Oat-mi the first one coming from solar neutrino analysis and the second one from atmospheric neutrino data. Since 
the solar electron anti-neutrino flux is strongly constrained by data [^fl I21I l^sl I2H fsoj l . that analysis puts severe limits 
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on the size of the magnetic moment (assuming a particular solar magnetic field profile) , in order to avoid strong spin- 
flavor conversion producing a sizable anti-electron neutrino flux. Recently also a limit for anti-non-electron neutrinos 
was quoted |it| . but these limits are weak and does not impose any constrain in our analysis. 

For the neutrino parameters in MSW-LMA region, the spin-flavor precession is very small for typical values of 
magnetic field in the Sun. However, it was recently pointed out |26j, |27j that random magnetic field could enhance 
this conversion. Consequently, stronger limits for the neutrino magnetic momentum /i was obtained, typically, 0.78 — 
1.2 x IO^Vb [3 

A conveniently chosen non- vanishing magnetic moment in the muon-tau sector leads to a very different scenario. Tau 
anti- neutrinos are produced through — > v r conversion, and assuming a vanishing mixing angle #13, the production 
of electron anti- neutrino is kept very small. The final solar neutrino flux can be a mixing of v e , u^, 2/ T , and u T , 
which can have some interesting phenomenological consequences. The more direct one would be a correlation between 
the solar magnetic field and the proportion between the different neutrino families. Also the regeneration effect will 
be modified due to a different proportion of active neutrinos in the solar mass eigenstates, in analogy with the effect 
of a non- vanishing #13 31] . The next round of reactor and long-baseline experiments can measure the #13 angle if 
sin 2 (26»i 3 ) > 0.01. 

We analyze here the scenario where neutrinos interact with random solar magnetic fields trough a non-vanishing 
magnetic moment between anti-muon and anti-tau-neutrinos as a sub-leading effect in the context of LMA solution to 
the solar neutrino anomaly. We combine the results of this analysis with the constrains coming from the KamLAND 
observations. 

II. FORMALISM 

We start working in a 6x6 matrix formalism, with v = (y e , v^, v T , v e , v^, v T ) T where we include, besides the usual 
mass induced oscillation, magnetic moment terms between second and third families. We use as the mixing matrix 
the standard PMNS (Pontecorvo, Maki, Nakata,Sakata) mixing matrix as presented in the Particle Data Group 
(PDG) [3]. After rotating out the angle 9 a t m = #23, we can decouple the first, second and sixth families, obtaining: 



-5c 2e + V e + V„ 5s 2 e 

{ ~d\ ( V Ji ) = ( 5s2e Scw + Vy. iiBexp(ia) | | < | , (1) 

^£?exp(— ia) A — 



where S = A ^ 21 , A = Am 32+^ Am 3i , Am| is the mass squared difference between neutrino families i and j, V e and 
are the matter potentials, a is a phase of magnetic field, C29 and S28 are cosine and sine of solar angle 9q. The 
eigenstates v' and v' T are linear combinations of weak states as v' = cg 23 ^ M + sg 23 v T , v' T = —Sd^v^ + C0 2a v T . (From 
now on we suppress the prime symbol.) 

It is more convenient to work in matrix density formalism, where the effects of random magnetic fields can be 
included in the evolution equation. The evolution equation, as given in Eq. can be rewritten using the formalism 



of the density matrix p, with elements pij = \vi >< = 1, ..,3. We can expand the resulting 3x3 matrix 

using a complete set of 3 x 3 matrices, : A^, {y = 0, 8): the Ao = \j2j 3/3 ( I3 is the 3x3 identity matrix) and 
\u(v — 1, ...,8) are the Gell-Mann matrices. We assume T^AyA^) = 2<5„ M . The final equation can be written as 



dt 



li, v, a = 0, . 



(2) 



where the h v = Trace(7?A 1/ )/2 are defined as H = Yl v =o 8 hv^vi where the elements h 2 ,h/i and /15 vanish. Similarly 
p = p v \ v . Explicitly the coefficients hi are 
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If Cav is identically zero, we have a Liouville equation, that after some algebra can be written as 
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where the matrix is antisymmetric as it is expected in the Liouville equation. 

Until now we are only dealing with the usual MSW mechanism with spin-flip terms via regular magnetic fields (see 
Eq. (IJ). Note that or he oc h 7 oc jj, — ► we recover the usual MSW LMA mechanism. 

A. Magnetic Field Profile 



In order to quantitatively perform the analysis, one has to choose a solar magnetic field profile. We assume for the 
magnetic field a triangular profile in the convective zone, with a maximum at r / ' Rsun = 0.85 of Bmax — 100 kG, and 
zero in the radiative zone. We assume that the magnetic field will be composed by a regular part and a random part. 
For values of oscillation parameters in the LMA region, the contribution of the regular field is completely irrelevant. 
However, the random character of magnetic field allows for large amplitudes changes that can significantly modify 
the neutrino evolution inside sun. 

We will introduce the random features in the magnetic field through a delta-correlated fluctuations (3^ , which 
has the advantage of allowing a simple analytical parametrization of the random field in the neutrino evolution 
equation. 



The system evolution can than be divided in two parts: first the simple MSW conversion in the production region 
of solar neutrinos, where we can take the formulas for the two families conversion, as presented, for instance, in |34j. 
For r I ' Rsun > 0.7 the magnetic field starts to act on the system and then the conversion probabilities will depend of 
the neutrino magnetic momentum fj,. 

The random features will be introduced through the piece of neutrino evolution in the matrix density formalism, 
Eq. J2J). In more general formalism, we should consider the relative size of coherent length of the magnetic field that 
we call Lq and the neutrino oscillation length, A„ = tt a a ■ The condition to have a decoupling between the LMA 
MSW oscillations and the spin-flip induced by random magnetic fields is defined as 

K » L Q . (5) 

Following Ref ^33] , after some algebra, the change in the neutrino evolution to the random character can be written 
as follows: 

Pll = P22 = P33 = Pii = P55 = /W 2 = P88/3 = ~2k, 

P38 = PS3 = 2V3fc , (6) 

k =< (pB x ) 2 > L =< (fJ-B y ) 2 > L , 

where B xy are the random components of the magnetic field perpendicular to neutrino trajetory. We take this 
components to be proportional to the regular magnetic field. 

To analyze if the values used for the parameter k are reasonable, we can write it in convenient units: 



where: 
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If we take the neutrino oscillation parameters from the best fit point of the standard solar neutrino analysis, 
(tan 2 9 Q ,Am%) = (0.4, 8 x 10~ 5 eV 2 ), we have that the oscillation length for a 10 MeV neutrino is A„ ~ 200 km. 

The probabilities can be treated classically since the averaging over the production region suppresses any interference 
effect. We calculate the final survival probability using: 
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(7) 



where Pli d is the probability that an electronic neutrino v e arrives at the bottom of the convective zone as v\ , P^° nv 
is the probability that a v\ crosses the convective zone and leave the sun as v e . Since we have B = in the radiative 
zone, P^ d = 0. The probabilities P™ d and P^ are equal to computed using the 2x2 evolution equation for the 
LMA MSW mechanism. 
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FIG. 1: Neutrino survival probabilities, the solid line is P ee , the dashed is the sum P ee + P efl , and the P e f is the remaining 
until the maximum. We used tan 2 8 = 0.4 in all panels. 

Now we are in condition to compute the Pf"™ and P^° nv probabilities. In the convective region, the regular 
magnetic field is too small to induce spin-flip conversion and then the elements h§ and /17 vanishes. Then the Eq. Q 
decouples in a sub-sector containing only (pi, p2j P3> Ps) elements as follows: 
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(8) 



It is interesting to note that the evolution equation does not depend on /i 8 , not depending therefore on the 
atmospheric mass scale Am| 2 - This is an important feature of our configuration, since as mentioned in the last sec- 
tion, the condition of validity of our equation is that the coherence length is smaller then the neutrino oscillation length. 



At the bottom of convective zone the initial conditions are 

v x : p {0) = 1/3, pi(0) = -sin2#/2, p 3 (0) = cos20/2, p 8 (0) = 1/(2^3) 
v 2 : p (0) = 1/3, Pl (0) = sin2#/2, p 3 (0) = -cos20/2, p 8 (0) = 1/(2^3) 



(9) 
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To study the anti-neutrino production, we solve numerically Eq. ||§J for different values of k, and with initial 
conditions given by the usual MSW mechanism in the radiative zone. 

The effect of the random magnetic field inclusion in evolution equation can be seen in Fig. ^ where we can read 
the fraction of the different neutrino flavors, v e : : v r . The solid line represents the electronic neutrino survival 
probability, P ee , while the dashed line is the v e + fraction, P ee + P efl The remaining until the no-oscillation value 
is the contribution of the P ef probability. In all panels we used tan 2 9 = 0.4. 

For small values of fc, the electronic and muonic neutrinos start to convert into v Tl and, as a consequence, the 
electronic survival probability decreases. For large values of k, the neutrino flux tends to split equally in the three 
neutrinos flavors, with 1/3 of the total flux for each flavor. 

Taking the marginal allowed value of Lo = 10 km, in order to have k = 10~ 14 eV as in the last panel of Fig. ^ 
we should have B ~ 10 MG, which is hardly acceptable. Lower values of neutrino energy will decrease the neutrino 
oscillation length, making it more difficult to fulfill the conditions in Eq. (JSJ. In this sense, the values of k = 10 -15 
eV seems more feasible in a realistic scenario. 

However, as pointed out in Section^ we expect that the effects of neutrino conversion would be stronger when A„ ~ 
Lq. So the limitations expressed in the last paragraph are only numerical limitations, and not physical constraints. 
Having this in mind, we decided to extend the analysis up to k = 10~ 14 eV, assuming that a complete numerical 
integration of Eqs. J3J and © would give qualitatively the same results as we present here. 

III. LMA REGION 

The validity of our numerical treatment is energy dependent, so we have made the choice to limit our fit to solar 
neutrino data to a specific neutrino energy range. Since the validity of our approximations may not hold for low energy 
neutrinos, we included only the data for the high energy neutrino (SNO-I, SNO-II and SK). We must be careful in 
analyzing allowed regions for neutrino parameters, since the inclusion of low energy solar neutrino data should change 
this picture. 

The effect of the random magnetic field in the LMA region can be seen in Fig. [3 As the parameter k increases, the 
electronic survival probability decreases. This effect can be compensated by a higher value of neutrino mixing angle, 
moving the allowed region to the right in left panel of Fig. [3 

For larger values of k, all probabilities tend to 1/3, with a weak dependence of the mixing angle. Since the 
Supcr-Kamiokande and SNO results are in accordance with this probability, in this scenario even maximal mixing is 
allowed. But also a interesting phenomenon occurs. Now a significant part of the total neutrino flux does not take 
part in regeneration effect in Earth. Actually, if we have exactly a equally equipartition of v ei and v T fluxes the 
regeneration effect vanishes, regardless of the neutrino mass difference Am 2 . 

The right panel of Fig. [2] presents the KamLAND allowed regions, for 95% C.L., 99% C.L. and 3<r. Maximal mixing 
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FIG. 2: The LMA compatibility region for SNO+SK data (left panel) and the KamLAND allowed region (right panel). In left 
panel, the black line stands for no magnetic field, the long-dashed line for k — 1CP 15 ' 5 , the short-dashed line for k = 1CP 15 , the 
dotted line for k — 10~ 14 ' 5 and the dot-dashed line for k — 10~ 14 . In the right panel the allowed regions stands for 95% C.L., 
99% C.L. and 3a for respectively the filled, dashed and black curves. 

is allowed at 62.1% C.L., and low values of Am 2 are still consistent with data at 99% C.L.. This last region is 
inconsistent with solar neutrino experiments because it predicts a too strong regeneration effect, not seen by data. 

In the present context, the regeneration effect can be suppressed for large values of k. As a result, lower values of 
Am 2 are allowed, and a new region of compatibility between solar and KamLAND data appears. In other contexts 
of non-standard neutrino physics [l^.[l^|. we have called this region very- low LMA. 

In Fig. |3 we present the allowed region for a combined analysis of high-energy solar neutrino and KamLAND data. 
We can see in this plot both the displacement of the allowed region to higher values of tan 2 9 for moderate values of 
fc, and the appearance of the very- low LMA region at Am 2 ~ [1 — 2] x 10 -5 eV 2 for k = 10 -14 eV. 

IV. 0i3 / 

In this work, we assumed #13 = 0, since non- vanishing values of this angle will lead to a production of 9 e , which is 
strongly constrained by data. Then, if #13 ^ a limit in /123 could be established in the same way the limits of the 



8 



10" 



> 

CD 



C\J 



10 



-5 




0.1 



1.0 



tg 2 e 



FIG. 3: Combined result for SNO+SK and KamLAND data, at 99% C.L., for the same values of k presented in Fig. [21 Maximal 
mixing is allowed for k > 1CP 14 ' 5 . 



other components of magnetic moments were found [26j], in a two families analysis, denoted here by \j?J am . 

Including the angle #13 in the mixing matrix would lead to the term [i v sin #13 connecting D e to the active neutrino 
families. So the limits in /j,?/ am could be scaled to a limit in /123 if a positive measurement of #13 is achieved in the 
next round of reactor |.3 r | and accelerator [3^ experiments. This limit would be weaker then the present limits in 



li j am by a factor 1/ sin( 



'13- 



CONCLUSIONS 



We have investigated new effects in solar neutrino phenomenology due to interactions of these particles with a 
random solar magnetic field. Considering the context of LMA realization of the MSW solution to the solar neutrino 
anomaly we have analysed the neutrino spin-flavor conversion phenomenon which appears as a subleading effect when 
a non-vanishing neutrino magnetic moment linking the second and the third families is assumed. Such a magnetic 
moment induces a large conversion of solar neutrinos into non-electron anti-neutrino flux, which is not severely 
constrained by the solar neutrino observations. 



9 

Since the mixing angle 6*13 can be considered very small, this conversion is not followed by a production of electron 
anti-neutrinos flux which, in contrast to non-electron anti-neutrinos, is very constrained by data. 

The results of our analyses of SNO+SK compatibility region indicate that in the presence of solar random magnetic 
fields the allowed region for Am 2 becomes larger while higher values of O12 arc found. 

This is a consequence of the fact that, in a three neutrino family context, an electron-neutrino survival probability 
of P ~ 1/3 is possible, even for tan 2 #12 ~ 1. In fact, when the random components of the solar magnetic field is large 
enough in such way that k = 10~ 14 , values around Am 2 w 10 _5 eV 2 are included in the region allowed by solar neutrino 
observations. Furtcrmore, different proportion of neutrinos and anti-neutrinos suppresses the regeneration effect of 
the solar neutrinos crossing the Earth matter. As a consequence, a totally new region of compatibility between solar 
neutrinos and KamLAND, which we call very-low LMA, appears at 99% C.L. for small values of Am^ <~ [1 — 2] x 10~ 5 
eV 2 and maximal mixing. 
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